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INTRODUCTION 

The  cylinder  is  a  shape  quite  frequently  encountered  In  engineering 
work,  and  the  nuclear  field  is  no  exception.   '  >ny  reactors,  many  fuel 
elements,  and  many  other  sources  of  radiation  have  shapes  approximating 
cylinders . 

It  is  not  uncommon  to  need  to  know  the  flux  from  the  end  of  such  a 
cylinder,  whether  the  flux  is  composed  of  neutrons  or  gamma  rays  or  both. 
This  information  is  required  for  such  reasons  as  the  design  of  adequate 
shielding  and  for  calculation  of  dose  rates.  To  the  author's  knowledge, 
there  is  no  method  of  calculation  available  at  this  time  which  can  be 
performed  without  a  computer  in  reasonable  time  periods  and  give  answers 
with  predictable  accuracy. 

An  approximate  solution  has  been  developed  by  Rockwell  (6)  for 
determining  upper  and  lower  limits  on  the  flux  from  the  end  of  a  cylinder. 
This  technique,  while  bracketing  the  flux,  in  many  cases  provides  such  a 
large  bracket  that  the  actual  flux  is  still  very  much  in  doubt. 

An  IEM  704.  Computer  code  has  been  developed  by  Gillis  et  al.  (3) , 
which  calculates  the  uncollided  flux  at  any  point  outside  a  cylinder. 

It  is  the  purpose  of  this  paper  to  develop  empirical  techniques 
enabling  one  to  determine  the  centerline  uncollided  flux  from  the  end  of 
a  finite  cylinder  without  the  aid  of  a  computer. 


NOMENCLATURE 

— P     1 
9  Scalar  flux  (cm-  sec"  ) 

— P    1 
S.        Source  strength  of  plane  source  (cm  sec"  ) 

Sy  Source  strength  of  volume  source  (cm  sec-  ) 

kj^  Macroscopic  cross  section  of  ith  shield  material  (cm  ) 

U  Macroscopic  cross  section  of  source  material  (cm-  ) 

t.  Thickness  of  ith  shield  material  (cm) 

tQ  Thickness  of  source  disk  (cm) 

a  Distance  from  end  of  cylinder  to  observation  point  (cm) 

_        n 

\i  2  p.t./a   (cm"1) 

i 

bi  J  Vi 

1 
z        Effective  self-attenuation  distance  (cm) 

b2        bl  +  ^s2 

h        Cylinder  height  (cm) 

b3  bx  +  Ush 

h'        Cone  height  having  same  volume  as  cylinder  (cm) 

y        *i +  ^h' 

R        Cylinder  radius  (cm) 

91  tan-1  RQ/a 

92  tan"1  RQ/(a  +  h) 

NMAX+1  Number  of  disks  in  cylinder 

6  Depression  of  equivalent  circular  plane  source  (cm) 

6N  tan"1  R0/(a+NtQ+6) 


62  tan"1  R0/(a+  z) 

E  (b)  b"'1  |     £l_  dt       n     >  0    ,     E  (b)  =  1— 


THEORETICAL  DEVELOPMENT 


Problem 


The  equations  for  calculation  of  uncollided  flux  from  some  geo- 
metrically simple  sources  are  developed  in  Rockwell  and  those  used  in 
this  work  are  listed  here.  From  a  circular  plane  source  (Fig.  l)  of 
radius  R  emitting  S.  particles  per  square  centimeter  -  second, 

0  (a)  =fA_JE1(b1)  -  E1(b1  sec  9)}  .  (l) 

From  a  truncated  cone  (Fig.  2)  of  height  h  emitting  Sv  particles  per 

cubic  centimeter- second  the  flux  at  P,  the  imaginary  apex  Of  the 

cone  is, 

SV     r  Mbi    sec©  )  E0(bo  sec  6)    ,      , 

*  <*> =  -  {Vbl>  "   J  9     -  Vb3>  +  2  3  fl    }  W 

«ip_  sec  o 

s 

The  flux  from  the  end  of  a  cylinder  (Fig.  3)  is  not  so  easily  ob- 
tained and  must  be  calculated  by  approximation  techniques.  There  are 
several  approximate  solutions  available  and  these  are  listed  and  dis- 
cussed. In  all  following  discussion  the  orientation  of  the  cylinder  will 
be  assumed  to  be  such  that  the  centerline  is  vertical  and  the  detection 
point  P  is  above  the  cylinder.  All  the  solutions  require  that  the  point 
P  be  on  the  centerline  of  the  cylinder.  All  solutions  also  require  that 
the  source  S^  be  of  constant  strength  throughout  the  cylinder. 

There  are  six  parameters  which  must  be  considered  in  this  problem 
area.  Those  associated  with  the  cylinder  are  Mg,  R  ,  and  h.  The 
parameters  associated  with  the  shield  are  b-,,  a,  and  _U  .  The  shield  is 
assumed  to  be  composed  of  slabs  perpendicular  to  the  cylinder  centerline. 


Fig. I.   Circular    plane   source.  Fig.  2.    Truncated    right- 

circular  cone  source. 
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Fig. 3.    Right -circular  cylindrical    source    viewed  endwise 


Solutions 

Stacked  Disk  Solution.  In  order  to  determine  empirical  relation- 
ships for  the  flux,  it  was  necessary  to  determine  the  flux  accurately 
and  this  solution  served  that  purpose.  A  cylinder  is  assumed  to  be  made 
up  of  a  series  of  thin  stacked  disks.  The  flux  at  P  from  each  disk  will 
be  calculated  by  assuming  each  disk  to  be  an  equivalent  circular  plane 
source  with  S  =  Syt  .  The  shielding  of  all  disks  above  the  one  under 
consideration  must  be  considered  so  each  disk  will  "see"  a  different 
amount  of  shielding. 

The  self-shielding  of  each  disk  is  considered  by  locating  the  plane 
source  some  distance  6  beneath  the  upper  surface  of  the  disk  as  shown 
in  Fig.  3.  Since  the  location  of  the  plane  source  within  the  disk  is 
actually  a  weighting  function  applied  to  the  source  distributed  within 
the  disk,  6  is  expected  to  be  less  than  t/2  because  the  upper  portion 
of  the  disk  is  more  important  due  to  the  fact  that  it  "sees"  less 
shielding. 

A  second  reason  exists  for  choosing  6  <  \>J2\   that  being  the  desire 

to  always  make  conservative  calculations  where  possible  exposure  to 

nuclear  radiation  may  be  involved.  As  the  plane  source  is  moved  closer  to 

the  upper  surface  of  the  disk,  it  "sees"  less  of  the  shielding  in  the  disk 

than  is  actually  the  case  so  the  resultant  flux  calculation  will  be  too 

high.  As  the  disks  are  made  thinner,  the  neglected  shielding  in  the  source 

will  decrease  and  the  solution  becomes  more  accurate. 

The  mathematical  expression  for  the  flux  is 
NMAX 

*  (a)=  2   N = '  0  {El(t>1  +Imst°  +  Ms6)  "  ^l+^o+^s)^0^}*  0) 


where  NMAX,  tQ,  and  6  will  be  chosen  to  assure  the  calculated  flux  is 
within  1%   higher  than  the  true  flux.  Equation  3  is  identical  to  Eq.  A-l 
in  Appendix  A  which  explains  the  computer  code  used  to  solve  the  equation. 

Large  Cone  Approximation.  This  solution  assumes  the  cylinder  to 
be  a  truncated  cone  the  same  height  as  the  cylinder  with  P  at  the 
apex.  The  angle  subtended  by  the  cone  is  9  on  Fig.  3,   and  the  cone  is 
obviously  larger  than  the  cylinder  it  represents.  The  flux  from  this 
approximation  is  calculated  by 


2U  L  2         sec  i       2      sec  e-. 


(4) 


This  solution  provides  an  upper  limit  on  the  flux  from  the  cylinder. 

Equation  A   is  identical  to  Eq.  A-3  in  Appendix  A  which  explains  the 

computer  code  used  to  solve  the  equation. 

Small  Cone  Approximation.  This  solution  also  assumes  the  cylinder 

to  be  a  truncated  cone  which  is  the  same  height  as  the  cylinder  with  P  at 

the  apex.  The  angle  subtended  by  the  cone  is  0p  on  Fig.  3,  and  the  cone 

is  obviously  smaller  than  the  cylinder  it  represents.  The  flux  from 

this  approximation  is  calculated  by 

c.         f                    E0(b-r  secOp)                        Eo(bo  sec  9J     -, 
*(a)=!L    {MbO     -S3 31    -E(b3)+2ll L    1        .      (5) 

2M   -  sec  9„      z ■        sec  0O    J 

s  <■  <c 

This  Solution  provides  a  lower  limit  on  the  flux  from  the  cylinder  and 
when  coupled  with  the  large  cone  approximation  provides  a  bracket  on  the 
expected  flux.  Equation  5  is  identical  to  Eq.  A-5  in  Appendix  A  which 
explains  the  computer  code  used  to  solve  the  equation. 

Equivalent  Volume  Cone  Appro ximat ion .  This  solution  also  assumes  the 
cylinder  to  be  a  truncated  cone  with  P  at  the  apex.  The  angle  subtended 


by  the  cone  is  G  on  Fig.  3,  and  the  height  of  the  truncated  cone  is  such 

that  the  volume  is  equal  to  that  of  the  cylinder  it  represents.  The  flux 

from  this  approximation  is  calculated  by 

a         (  E_(b.   secO)  E„(b'   sec  8.)  1 

ft  (a)=  3L     JE  (bj   -  -i-i L.  -  E  (tf  )  *  2J i-  }     .       (6) 

2H  sec  8  sec  9, 

si  1 

Equation  6  is  identical  to  Eq.  A-4-  in  Appendix  A  which  explains  the 
computer  code  used  to  solve  the  equation. 

Equivalent  Circular  Plane  Source  Approximat ion .  This  solution 
assumes  that  the  cylinder  can  be  represented  by  a  single  circular  plane 
source  the  same  radius  as  the  cylinder.  The  equivalent  source  is 
assumed  to  be  located  somewhere  within  the  space  occupied  by  the  cylinder 
it  represents  so  that  the  shielding  it  "sees"  is  that  of  the  cylinder 
within  the  confines  of  the  cylinder.  To  allow  for  self-shielding  within 
the  source,  the  strength  of  the  equivalent  source  must  be  considered  in 
two  parts.  For  cylinders  with  small  mean  free  path  height  (|igh)  the 
equivalent  source  strength  will  be  assumed  to  be  S.  =  Syh.  For  cylinders 
with  a  mean  free  path  height  greater  than  some  value  Kj,  which  is  chosen 
empirically,  the  source  strength  will  be  assumed  to  be  S.  =  KSy  .  The  flux 
from  this  approximation  is  calculated  by 


"5 


V    r  i 

ft  (a)  =  JL    |  E1(b2)  -  jSj(b2  sec  9  J   |,   if  Mi  <K. 


(7) 


or 


KS. 


ft   (a)  =  — ?    {  E1(b2)  -  E1(b2  sec  0J    [,   if  Jigh>K.  (8) 


2^s 


Equations  7  and  8  are  identical  to  Eqs.  A-6  and  A-7  respectively  in 
Appendix  A  which  explains  the  computer  code  used  to  solve  the  equations. 


NUMERICAL  ANALYSIS 
Method  of  Solution 

This  problem  has  six  parameters  which  may  vary;  R  ,  U  ,  h,  U,  b, ,  and 
a.  In  order  to  consider  all  possible  combinations  of  parameters  adequately, 
a  large  number  of  problems  had  to  be  solved.  To  do  this  an  IBM  650 
Computer  code  was  written  which  solved  the  problems  by  each  of  the  five 
previously  listed  solutions.  This  code  is  described  in  Appendix  A. 

By  proper  choice  of  disk  thickness  the  stacked  disk  solution  could 
be  used  to  solve  the  problem  with  any  desired  degree  of  accuracy.  This 
solution  was  used  to  provide  an  answer  accurate  to  within  1%.     The  choice 
of  disk  thickness  to  assure  this  accuracy  will  be  discussed  later. 

The  three  conical  approximate  solutions  were  compared  with  the 
stacked  disk  solutions  and  the  relative  error  of  each  solution  was  cal- 
culated. The  results  were  examined  in  an  effort  to  determine  some 
simple  method  of  empirically  presenting  the  data  to  enable  one  to  calculate 
the  correct  flux  within  some  set  limits  of  accuracy  without  the  use  of 
a  computer. 

The  plane  source  approximation  was  used  as  a  trial  and  error  cal- 
culation to  determine  the  self-absorption  distance  which  forced  the  solution 
to  approximate  the  stacked  disk  solution  within  specified  limits.  The  results 
were  examined  in  an  effort  to  determine  a  simple  empirical  relationship 
between  the  self-absorption  distance  and  the  flux. 

Preliminary  Investigations 

Accuracy  Analysis  of  Stacked  Disk  Solution.  Using  the  computer  code, 
it  was  first  necessary  to  determine  the  disk  thickness  which  gave  sufficient 


accuracy.  As  the  disks  became  thinner  the  solution  became  more  accurate, 
so  very  thin  disks  were  preferred.  On  the  other  hand,  computer  running 
time  increased  as  disk  thickness  decreased  since  more  disks  were  then 
required  for  any  given  problem.  A  compromise  was  made  such  that  the  flux 
was  always  calculated  to  within  1%   of  the  asymptotic  value,  the  result 
being  always  higher  than  the  true  value.  The  required  disk  thickness 
could  then  be  determined  for  different  combinations  of  parameters. 

In  order  to  better  understand  the  problems  associated  with  the  work 
undertaken,  it  was  necessary  to  approach  the  problems  of  disk  thickness 
and  location  of  the  equivalent  plane  source  within  each  disk  in  an 
oblique  fashion.  It  was  first  assumed  that  the  plane  source  was  on  the 
upper  surface  of  each  disk,  i.e.  6=0,  and  disk  thicknesses  were  deter- 
mined. Problems  solved  with  this  choice  of  6  were  quite  lengthy,  due  to 
the  necessity  of  using  extremely  thin  disks  since  no  self -absorption  in 
each  disk  was  considered.  This  extreme  location  of  6  was  used  initially 
since  the  results  would  satisfy  the  original  requirement  that  the  flux 
always  approach  the  correct  flux  from  some  higher  value. 

With  the  information  provided  from  the  initial  set  of  problems,  it 
was  possible  to  compare  the  results  with  those  obtained  when  6  was  varied. 
Three  facts  were  immediately  apparent  when  these  comparisons  were  made. 
First,  the  optimum  location  for  6  was  quite  near  t  /2.  Second,  the 
optimum  location  for  6  varied  as  the  parameters  of  the  cylinder  and 
shield  varied.  Third,  the  major  portion  of  the  advantage  of  the  depression 
of  the  plane  source  was  observed  as  6  varied  from  0  to  0.4.t  .  The  choice 
was  made  for  6=  0.4-5tn  and  proved  to  be  satisfactory.  The  advantage  of 
this  depression  of  the  plane  source  was  a  decrease  in  the  number  of  disks 
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required  for  any  given  cylinder,  the  decrease  being  at  least  a  factor  of 

ten  which  saved  considerable  computer  time. 

In  the  search  to  determine  the  proper  disk  thickness  for  each  problem, 

three  parameters  were  shown  to  have  an  appreciable  effect  on  the  thickness 

necessary  to  provide  the  1%   accuracy  required.  The  results  of  this  search 

are  shown  in  Table  1  with  the  variation  of  a  single  parameter  from  some 

intermediate  value  shown  as  the  cause  with  the  change  of  disk  thickness 

while  still  maintaining  1%   accuracy  of  the  calculated  flux  shown  as  the 

effect . 

Table  1.  Effect  on  disk  thickness  required  to  determine  flux 
within  1%  of  the  asymptotic  value  as  parameters  are 
varied  from  intermediate  values. 

Cause  :  Effect 

■ 

Increase  b-i  Thicker 

Increase  RQ  Thicker 

Increase  M-s  Thicker 

The  disk  thickness  was  determined  in  terms  of  mean  free  paths,. 
Mst  .  Use  of  mean  free  paths  makes  it  easier  to  explain  the  dependence 
on  JJg  while  not  affecting  the  other  results. 

By  the  choice  of  6,  the  location  of  the  circular  plane  source  within 
each  disk  representing  the  disk,  some  small  amount  of  the  shielding  within 
each  disk  has  not  been  considered.  This  was  done,  as  previously  explained, 
to  force  the  flux  calculation  to  always  approach  the  asymptotic  value  from 
some  higher  value.  Certainly  one  expects  this  neglected  shielding  to 
become  less  important  as  the  total  amount  of  shielding  is  increased.  This 
expectation  is  borne  out  by  the  observed  effect  that  increasing  b-,  has  on 
disk  thickness. 
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In  considering  the  dependence  on  RQ>  it  is  only  logical  to  assume  that 
to  represent  a  disk  by  a  plane  the  disk  must  "look  like"  a  plane,  that  is, 
RQ  »  tQ.  This  being  the  case,  one  expects  that  disks  with  larger  radii 
could  be  thicker  and  still  "look  like"  a  plane.  This  expectation  is  boruu 
out  by  the  observed  effect  that  increasing  R  has  on  disk  thickness. 

The  effect  of  Hs  on  disk  thickness  is  again  one  of  geometry.  As  p. 
increases,  the  disk  thickness  in  terms  of  mean  free  paths  is  expected  to 
increase  because  this  tends  to  hold  the  shape  of  the  disk  constant.  If  the 
thickness  did  not  increase,  the  disk  would  approximate  a  plane  even  more 
closely  than  necessary.  The  expected  results  were  observed  and  disk 
thickness  did  increase  as  Pg  increased. 

Parameter  Analysis.  Having  chosen  6  it  was  possible  to  calculate  the 
flux  from  cylinders  whose  parameters  were  chosen  at  random,  as  were  the 
parameters  of  the  shields.  These  problems  were  solved  in  an  effort  to 
determine  which  approximate  method  of  solution  or  combination  of  methods 
showed  the  greatest  promise  as  a  tool  to  be  eventually  used  in  the  simple 
empirical  presentation  of  data  which  was  desired. 

The  three  conical  solutions  were  shown  to  have  little  promise.  As 
^s^o  decreased  no  single  solution  or  average  combination  of  the  approximate 
conical  solutions  could  provide  accurate  results.  Errors  were  greater 
than  a  factor  of  10  in  some  cases.  The  errors  varied  in  such  inconsistent 
fashion  that  no  pattern  which  might  make  the  application  of  correction 
factors  to  the  approximate  solution (s)  feasible  could  be  observed.  In 
most  cases  the  large  cone  approximation  and  the  equivalent  volume  cone 
approximation  were  more  nearly  correct  and  the  large  cone  approximation  was 
the  more  reliable  of  these  two  because  it  always  gave  a  conservative  result. 
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The  equivalent  plane  source  approximation  was  chosen  as  the  only 
approach  which  showed  some  promise  of  lending  itself  to  the  empirical  fit 
of  large  amounts  of  data.  It  was  observed  that  for  large  |isR  ,  the  self- 
absorption  distance  in  terms  of  mean  free  paths,  UgZ,  depended  only  on 
bn  for  a  given  cylinder  height  and  it  was  decided  to  present  the  results 
as  a  family  of  curves  for  different  ugh  with  ugz  plotted  versus  b^. 

It  was  possible  at  this  point  to  choose  a  value  for  K,  the  mean  free 
path  cylinder  height  that  fixed  the  constant  source  strength,  for  cylinders 
whose  height  exceeded  K  mean  free  paths.  The  artificial  circular  plane 
source  having  S.=  Sy  h  kept  increasing  as  h  increased,  although  the  flux 
from  the  cylinder  did  not.  The  self-absorption  distance  z  necessarily 
increased  as  the  plane  source  increased  and  these  compensating  effects 
were  undesirable  since  the  flux  did  not  continue  to  increases.  The  value 
of  K  was  fixed  at  1.5  so  that  all  cylinders  more  than  1.5  mean  free  paths 
tall  were  assumed  to  have  an  equivalent  circular  plane  source  strength  such 
that  S.  =  1.5  Sy/|i .  The  plane  source  strength  then  fell  into  one  of  two 
categories:   if  (I  h  <  1.5,  S^  =  Svh;  if  p.  h  >  1.5,  SA  =  1.5  Sy/us. 
The  figure  of  u  h  =1.5  was  chosen  as  the  minimum  value  which  could  be 
used  to  fix  a  constant  source  strength  without  a  risk  of  developing  a 
negative  z,  which  is  not  anticipated  in  the  computer  program.  The  minimum 
value  was  desired  since  the  accuracy  of  the  empirical  relationships  improved 
as  the  cylinder  height  used  to  fix  a  constant  source  strength  decreased. 

Data  Analysis 

Typical  values  for  RQ  and  Ug  were  chosen  and  problems  with  varying 
±f>  .§»  bl>  a-11^  h   were  solved  and  the  resultant  curves  in  Figs.  4a  and  4b  were 
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determined.  The  self-absorption  distance  Ugz,  expressed  in  terms  of  mean 
free  paths,  refers  to  the  equivalent  circular  plane  source  which  is 
assumed  to  replace  the  cylinder.  Presentation  of  the  data  in  this  manner 
requires  that  Ugz  be  dependent  only  on  the  shield  mean  free  paths,  b^, 
for  a  given  cylinder  height.  This  requirement  is  satisfied  for  large 

For  the  cases  of  small  U  R  ,  marked  deviation  of  URz  from  that  sug- 
gested by  Figs.  4-a  and  4-b  was  observed.  In  general,  Ugz  increased  as 
UgR  first  dropped  below  some  critical  value  dependent  on  b-^,  and 
eventually  decreased  rapidly  as  U  R  continued  to  decrease.  Figure  5 
indicates  the  area  where  fluxes  may  be  calculated  to  within  10$  of  the 
actual  value  without  the  use  of  correction  factors.  As  (J.gR0  exceeds  the 
critical  value  shown  for  any  b^,  correction  factors  are  not  required. 

The  maximum  deviation  of  Ugz  from  the  values  predicted  by  Figs.  4-a 
and  4b  is  observed  for  large  jU.  For  a  given  value  of  UgR0  the  deviation 
was  greater  for  smaller  Ug  as  _U  was  held  constant.  The  relative  deviation 
was  slightly  dependent  on  U  h  and  quite  dependent  on  b^.  It  was  necessary 
to  develop  correction  factors  for  Jjlsz  which  depend  on  b-,,  JJ-,  and  Ug.  This 
has  been  done  in  Figs.  6  through  12.  These  correction  factors  were  deter- 
mined by  assuming  cylinders  having  U  h  =  2.  and  the  results  are  most  accurate 
for  this  case.  For  values  of  U  h  not  equal  to  two  the  error  in  the  cal- 
culated flux  for  small  UgR0  increases.  Largest  errors  are  noted  when  Ug 
is  small,  R  is  small,  and  u  is  large.  Errors  are  greatest  for  the  smallest 
Ush.  Errors  are  positive,  i.e.  the  empirically  determined  flux  is  greater 
than  the  flux  determined  by  the  stacked  disk  solution,  for  most  cases, 
and  the  estimated  maximum  negative  error  is  <25$.  The  estimated  maximum 
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positive  error  is  <  115%   and  these  extreme  errors  apply  only  in  isolated 
cases.  Table  2  lists  the  estimated  maximum  errors  anticipated  for 
varying  |ig  and  varying  Hsh  at  values  of  Hsh  other  than  2.  The  comparable 
errors  resulting  from  the  use  of  the  large  cone  approximation  are  included 
in  Table  2.  Appendix  B  explains  the  computer  program  used  to  calculate 
the  error  in  flux  determination  by  the  empirical  method. 
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DESIGN  CONSIDERATIONS 

Use  of  Design  Curves 

The  following  curves  are  to  be  used  to  determine  the  location  or 

self-absorption  distance  of  a  circular  plane  source  assumed  to  represent 

a  cylinder  for  purposes  of  calculating  the  centerline  uncollided  flux  at 

some  point  P,  at  distance  a  beyond  an  end  of  the  cylinder.  The  self- 

absorption  distance  is  measured  from  the  end  of  the  cylinder  closest  to 

P.  Necessary  input  data  are  ji,  a,  ps,  R,-,,  h,  and  b1.  The  following 

step-by-step  procedure  enables  one  to  determine  the  self-absorption  distance 

and  the  resultant  flux.  The  following  limits  on  the  values  of  the 

parameters  must  be  observed: 

R  >  1  cm 
0  - 

O.lSb-,^  <  A.0 
0.02cm"1  <  |J.3 
0.02cm   S  U 
To  determine  the  self-absorption  distance  proceed  as  follows: 
Step  1.  Calculate  Ugh. 

Determine  p.  z   from  Fig.  4a  or  4-b. 

Step  2.  Calculate  p.  R  . 

r  so 

From  Fig.  5  determine  necessity  of  correction  factor 
application.  If  correction  factor  is  necessary  go 
to  step  3;  if  not,  go  to  step  6. 
.  Step  3.  Using  Figs.  6b,  7b,  8b,  9b,  10b,  lib,  or  12b 

select  the  Figures  at  the  b,  values  immediately  larger 
and  smaller  than  the  given  b, .  From  these  two  Figures 
determine  a  relation  to  the  appropriate  F^,  G^,  %,  J^, 
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K.,  L. ,  or  M-  curves.  For  example,  at  b-,  -  U,   on 
Fig.  8b  at  H  =  0.3  and  u  -  0.1  a  location  between 
H?  and  H_  is  observed  approximately  one-third  of  the 
distance  from  H„  towards  H^. 

Step  U'     With  the  relations  determined  in  Step  3,  use  Figs  6a, 
7a,  8a,  9a,  10a,  11a,  or  12a  as  appropriate,  and  by- 
interpolation  determine  correction  factors  at  the  two 
values  for  b-^  used  in  Step  3. 

Step  5.  Interpolating  between  b-^  values,  determine  correction 
factor,  CF,  to  be  used  at  the  given  b-^,  and  multiply 
CF  by  U  z  which  is  then  the  M-Sz  to  be  used  in  succeeding 
steps. 

Step  6.  Calculate  sec  9Z  and  z  from  ugz, 


where  sec  6 


-M 


o  _ 
2 


(a+z) 
Step  7.  Calculate  b  =  b-,  +  Hsz. 

Calculate  \E1(b2)  -  E]_(b2  sec  8JJ  , 

Y2    3 

where  E^x)  =  -0.5772  -  In  x  +x  -  j-  +  |- 
if  x  <1,  Refs.  (1,2,6). 


(x)  =  e^  j  0.251+  2.335X*  x^  j? 
x    1.082+"3.331x  +  x'i 
if  x_>l,  Ref.(3). 


Step  8.  If  Hsh<1.5,  SA  =  Syh 

If  Msh>1.5,  SA  =  Sv(i^) 


^s 


Step  9.  «  (a)  =  -i  {  E1(b2)  -  E1(b2  sec  0z)j 


2 

If  the  flux  is  desired  at  some  point  P  not  on  the  centerline  of  the 
cylinder  a  conservative  estimate  can  be  obtained  by  assuming  the  point  P 
is  on  the  centerline  of  the  cylinder  at  the  same  vertical  distance  above 
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the  source,  and  solving  the  problem  in  the  manner  outlined  above. 

Sample  Problem 

The  following  problem  will  serve  to  illustrate  the  use  of  the 
design  curves.  The  parameter  values  have  been  assumed  to  be  as  follows: 
\is   =  0.07,  ji  =  0.3,  hi   =  6,  RQ  =  10  cm,  h  =  4.0  cm,  a  =  20  cm. 
Step  1.  Hsh  =2.8 

From  Fig.  4-b,  UgZ=  0.497 
Step  2.  HgR0  =0.7  From  Fig.  5,   correction  factor  is  necessary. 
Step  3.  From  Fig.  8b,  observe  a  location  between 

H_  and  Ho  close  to  H.,. 

From  Fig  9b,  observe  a  location  between 

J*  and  Jp  close  to  Jp» 
Step  4-.  From  fig.  8a,  by  interpolation,  the  correction  factor 

is  estimated  as  0.97.  From  Fig.  9a,  by  interpolation, 

the  correction  factor  is  estimated  as  1.01. 
Step  5.  Correction  factor  is  0.983  for  b^^  =  6. 

Corrected  kz  is  0.4-89. 
Step  6.  z  =  6.96  cm. 

sec  0„  =  1.0665. 
Step  7.  b2  =  6.4-89. 

b2sec  Qz   =  6.920. 

E1(6.4-89)  =  2.078  x  10"4". 

E1(6.92l)  =  1.273  x  10"^. 
Step  8.  SA  =  21.4.28  cm-2  sec 
Step  9.  $(a)  =  8.625  x  10-^  cm"2  sec"1 
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1-2         -1 
The  stacked  disk  solution  to  this  problem  is  7.8  x  10"^  cm   sec  , 

The  error  in  the  empirical  solution  is  10. 656.  The  error  in  the  large 

cone  approximate  solution  is  64»7$. 
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DATA  PRESENTATION 
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Fig.  6b.    Source    and    shield     cross     sections       correlation    for   b,=  0.1. 
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Fig.  7b.    Source      and     shield      cross     sections       correlation    for    b,=  I.O. 
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Fig.  8b.    Source      and     shield      cross     sections      correlation    for    b,»  4.0. 
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Fig.  9b.    Source     and    shield     cross     sections       correlation    for   b,=  10. 
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/X8     SOURCE      CROSS      SECTION    (cm.-') 


Fig.  |0b.    Source     and    shield     cross     sections       correlation    for    b,=  20 
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jJL%    SOURCE      CROSS      SECTION    (cm:') 


Fig.  Mb.    Source      and     shield      cross     sections      correlation    for    b,=  30. 
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Fig.  12b.    Source     and    shield     cross     sections       correlation    for    b,=  40. 
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CONCLUSIONS 

The  largest  errors  in  determination  of  fluxes  by  use  of  the  curves 
presented  in  this  paper  will  occur  in  the  region  where  correction  factors 
are  necessary,  due  to  the  interpolation  required  to  determine  correction 
factors.  Even  so,  this  is  the  region  where  the  greatest  benefit  is 
derived  from  this  work,  for  it  is  here  that  currently  used  conical  approxi- 
mations are  most  in  error. 

For  any  given  cylinder  height,  h,  the  cylinder  ceases  to  look  like 
a  disk  and  becomes  a  slender  rod  as  R0  decreases.  Similarly,  for  a  fixed 
mean  free  path  height,  ^sh,  as  Jig  decreases  the  cylinder  becomes  longer 
and  for  a  fixed  radius  the  cylinder  again  appears  rod- like.  Certainly 
the  artificial  nature  of  the  plane  source  approximation  would  be  expected 
to  have  adverse  effects  in  this  region,  if  at  all,  and  the  correction 
factors  necessary  for  small  HgR0  were  no  surprise  although  their  shape 
could  not  be  anticipated. 

An  empirical  solution  similar  to  the  one  presented  herein  is  presented 
in  Rockwell  (6)  which  enables  one  to  determine  the  flux  from  the  side  of  a 
cylinder  by  approximating  the  cylinder  by  a  line  source  located  within  the 
cylinder.  This  empirical  solution  provides  a  method  amenable  to  hand  calcu- 
lation to  reproduce  the  results  shown  by  Taylor  and  Obenshain  (7) .  This 
original  work  is  applicable  only  for  cylinders  having  h»  R  and  h>  VVS» 
These  limitations  are  not  mentioned  by  Rockwell  and  therefore  his  empirical 
solution  may  be  frequently  applied  in  instances  which  exceed  the  restrictions 
placed  on  the  original  data. 

It  is  recommended  that  work  be  done  in  the  area  of  accurate  flux  solutions 
from  the  side  of  a  cylinder  for  cases  excluded  by  Taylor  and  Obenshain  (7), 
and  that  an  effort  be  made  to  combine  the  results.  Further  effort  should  be 
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made  to  present  these  results  in  some  empirical  form  which  lends  itself  to 
hand  calculation.  In  particular,  it  is  recommended  that  an  effort  be 
made  to  present  the  results  in  a  manner  analogous  to  that  used  in  this 
paper. 
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APPENDIX  A 


Description  and  Explanation  of  the  IBM  650  Computer  Program 
Used  to  Calculate  the  Centerline  Uncollided  Flux 
from  the  End  of  a  Finite  Right-circular  Cylinder 


This  program  was  written  to  calculate  the  centerline  uncollided  flux 
from  the  end  of  a  finite  cylinder,  both  accurately  and  by  existing 
approximation  methods,  and  to  calculate  the  self-absorption  distance  of 
an  equivalent  source.  This  program  was  written  in  Soap  II  computer 
language  using  floating  point  operations.  The  logic  diagram  and  the  object 
program  are  given  in  this  appendix. 

The  accurate  solution  is  designated  FLUX  and  the  solution  programmed 

was 

UMAX 


(a)  s  !_2      I     {E1(b1  +(NtO.«)  tyg  -  £^(^+(^0.45)^)   sec  6N  ]} 

2  f"°    \ U-i) 

where  eee  9n  =      J*  _____  2  ^ 


Equation  A-l  is  derived  as  Eq.  3,  the  stacked  disk  solution  on  page  5,  with 
6  set  at  0.45. 

The  solution  to  the  large  cone  approximation  was  designated  FLUXB 

and  the  solution  programmed  was 

S   f  Bl(K  sec  8.)  E(b  sec9}". 

•  «._.{  v*i> "  219       "  EW  +      Lee     ' •      C  3) 

2^.    ^         sec  o,  sec  ui 

Equation  A-3  appears  as  Eq.  4-  on  page  6. 


4o 


The  solution  to  the  equivalent  volume  cone  approximation  was  designated 

FLUXM  and  the  solution  programmed  was 

S„       r  E  (b     sec  9,)  E  (bi   sec  Q_), 

*  (a)   =  2-  {  E?(bl)   -  J-4 L  -  E2(b')+  ±1 11}  (A-4.) 

2US  2  sec  91  2     J  sec  G;L      J 

Equation  A-4.  appears  on  page  7  as  Eq.  6. 

The  solution  to  the  small  cone  approximation  was  designated  FLUXS 

and  the  solution  programmed  was 

S       r  E  (b     sec  G   )  E  (b.   sec  6  J     1 

<b  (a)  =  JL  {  E  (bj   -  J-l L  -  B.0O+  _2_2 2_     }   .       (A-5) 

2US         2     X  sec  G2  2     J  sec  G2  J 

Equation  A-5  appears  on  page  6  as  Eq.  5. 

The  solution  to  the  equivalent  circular  plane  source  approximation 

was  designated  FLUXZ  and  the  solutions  programmed  were 

<£  (a)  =  f3-  JE1(b2)  -  Ej  (b2  sec  9Z)  }  ,  if  Hsh  <  1.5  (A-6) 

and 

SV(1.5)  f  , 

<*>  (a)  =  -I \\(h2)   -  EiCbg  sec  6a)J  if  ugh>  1.5.  (A-7) 

2^s 
This  was  a  trial  and  error  solution  where  z  was  altered  as  necessary  to 

cause  FLUXZ  to  approximate  FLUX  within  a  specified  accuracy.  Equations  A-6  and 

A-7  appear  on  page  7  as  Eqs.  7  and  8  with  K  set  at  1.5. 

The  program  subrouting  for  E, (x)  used  different  solutions  for  values 

of  x  <1  and  for  x_>  1. 

For  x<l 

E,(x)  =  -0.57721566  -  In  x  +  x  -  x2. +x£  _  ...  .  (A_g) 

A      18 

For  x  >  1 

Mx)  =  £f  {  a°  1  ai  x  *  a2 1. 1  *3  x3  -  t. }  (A-9) 

x   bQ  +  bi  x  +  b2  x2  +  b3  x^  +  x4- 
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where 

ao 

=  0.26777373                 bQ  -   3.958^969 

al 

=  8.6347609                  \  =  21.099653 

a2 

=  18.059017                  b2  =  25.632956 

a3 

=  8.5733287                  b3  =  9.5733223 

Equatior 

1  A-8  is  an  infinite  series  and  the  choice  of  number  of  terms  used 

was  set 

by  requiring  that  the  ratio  of  the  absolute  value  of  the  last 

term  to 

the  complete  series  be  less  than  some  specified  value.  Equation  A-9 

is  an  approximate  solution  developed  by  Hastings  (5),  accurate  to  Within 

the  limitations  of  the  computer. 

Input  and  output  data  associated  with  this  program  are  listed  in 

Tables  A-l  and  A-2  respectively. 

Table  A-l.  Input  data  required  for  use  of  the  IBM  650 

Computer  program  which  calculates  the  centerline 

uncollided  flux  from  the  end  of  a  finite  right- 

circular  cylinder. 

Symbol 

Explanation                       Drum  Storage 

Location 

NMAX 

One  less  than  the  number  of  disks  in  the  cylinder      0100 

TZRO 

Disk  thickness  (cm)                              0101 

MUS 

Cylinder  cross  section  (cm-1)                      0102 

BONE 

Shield  mean  free  paths                            0103 
Shield  thickness  C  cm)                             0104 

A 

R 

Cylinder  radius  (cm)                             0105 
Source  strength  (cm-- '  sec"  )                       0106 

S 

CHK 

Check  to  calculate  or  pass  by  a  given  cylinder        0107 

ZZRO 

Initial  self -absorption  values  (cm)                  0108 

CRIT 

Required  precision  of  equivalent  source  flux          0500 

calculation 

CRIT1 

Cylinder  height  (mean  free  paths)  for  constant        0550 

equivalent  source  strength 

PT45 

Ratio  of  6  to  disk  thickness                       0350 
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Table  A-2.  Identification  table  for  print  outs  for  center line 
uncollided  flux  calculation  from  the  end  of  a 
finite  cylinder. 

Card  Word  Symbol  and/or  Explanation 

1  1  FLUX 

2  FLUXB 

3  FLUXM 
U  FLUXS 

5  FLUXZ 

6  z  Self-absorption  distance  (cm) 

7  (FLUX  -  FLUXZ)  /  FLUX 

8  M'gZ  Self-absorption  distance  mean  free  paths 

2  1  FLUXB  -  FLUX 

2  (FLUXB  -  FLUX) /FLUX 

3  FLUXM  -  FLUX 

U  (FLUXM  -  FLUX) /FLUX 

5  FLUX  -  FLUXS 

6  (FLUX  -  FLUXS)/FLUS 

3  1  HITE  Cylinder  height  (cm) 

2  NMAX 

3  TZRD 
U  MUS 

5  BONE 

6  A 

7  R 

8  S 

4  1  HITE+  A 

2  (HITE)  (MUS)  Cylinder  height  mean  free  paths 

3  (TZRD)  (MUS)  Disk  thickness  mean  free  paths 
4.  (R)  (MUS)  Radial  mean  free  paths 

5  MU  Shield  cross  section  (cm-1) 

6  A+R 

7  A/R 

8  A/HITE 


Approximate  running  time  for  the  program  is  5  minutes  for  NMAX  equal  to 
100.   The  program  is  faster  when  b-.  >  1  due  to  the  difference  in  the  series 
calculations  to  be  made. 


a.  tc 
<t  < 


u 

X  ■ 

O  (9 
2  <   < 

P  U  UJ 


<  z 
o  o 

<  ■ 


r> 


*  o 

-I 
5  o  o 


_i  _ 

<  „7 


IT 

<     „ 


_l        « 

<    m 


cm  y  z 
2;  < -1 


3      o 

O  -".CO 


~j«J  0|iJ 

a  tc  <» 


*.?K 


3  z 


z 

4 
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o  z 
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ll 
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HEIGHT  0 
ATED     COI 
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I 

o 

UJ 
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o 
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< 
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s 
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> 

o 

_j 
< 

TRUNC 
HAVINC 
VOLUM 
CYLIN 

' 

-      z 
o-  ►"  7 

O  111    I 

J»N 


o  m     z 
_,co     i 


I    = 


<  ?  "• 

S  i  i 


a  x 

°3 


X 

5  2 


H     o 


10 


Q.  Z  B) 
O  3  Ul 

Jilt 


A  U- 

H  I 
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O  X 
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o  —  a: 

<  <t  a. 

u  l  + 
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OBJECT  PROGRAM  FOR  CALCULATION  OF  FLUX  AND 
EQUIVALENT  CIRCULAR  PLANE  SOURCE  SELF-ABSORPTION 
DISTANCE 


BLR 

1951 

I960 

1 

0000 

00 

0000 

OOOO 

BLR 

1977 

19B4 

3 

0000 

00 

000  0 

OOOO 

BL  R 

1700 

1749 

3 

0000 

00 

000  0 

OOOO 

3  Y  N 

■  mi 

0100 

0I8K8LE881 

4 

0000 

00 

0000 

OOOO 

8  V  N 

T  Z  RO 

0101 

0  1  8K  TH  KNE  B 

5 

0000 

00 

0000 

OOOO 

8  Y  N 

MU  8 

0103 

SOURC  E   A  e 9 

6 

0000 

00 

0000 

oooo 

8  Y  N 

BONE 

0103 

SHIELD   WFR 

7 

0000 

00 

0000 

oooo 

8  Y  N 

A 

0104 

SHIELD   THK 

B 

0000 

00 

0000 

000  0 

S  Y  N 

R 

0105 

CYL   RADIUS 

9 

0000 

00 

000  0 

oooo 

8Y  N 

8 

0106 

8  0  II  R  C  E 

10 

0000 

00 

oooo 

oooo 

8  Y  N 

CH  K 

0107 

P08   OR   NEG 

11 

0000 

00 

0000 

oooo 

8YN 

ZZRO 

0108 

13 

0000 

00 

oooo 

oooo 

SYN 

START 

1000 

START   LOC  « 

13 

0000 

00 

oooo 

oooo 

ZK  0 

00 

0000 

0000 

ZER  0 

14 

0000 

00 

oooo 

oooo 

ONE 

10 

0000 

0051 

ONE 

15 

0050 

10 

oooo 

0051 

T*  0 

30 

0000 

0051 

TWO 

16 

0150 

30 

oooo 

0051 

THREE 

30 

0000 

0051 

THR  EE 

17 

0200 

30 

oooo 

0051 

FOUR 

40 

oooo 

0051 

FOUR 

IB 

0250 

40 

oooo 

0051 

TEN 

10 

0000 

0053 

TEN 

19 

0300 

10 

oooo 

0052 

RT45 

45 

0000 

00  50 

POINT   45 

30 

0350 

45 

oooo 

0050 

FP  100 

10 

0000 

0053 

HUNOR  E  D 

31 

0  400 

10 

oooo 

OOS  3 

FPBO 

S  0 

0000 

0053 

EIGHTY 

22 

04  50 

80 

oooo 

0053 

CR  1  T 

5  0 

0000 

0047 

Z   CRITERIA 

33 

0500 

50 

000  0 

0047 

CR  1  Tl 

IS 

0000 

0051 

24 

0550 

15 

oooo 

0051 

P  1 

31 

4  159 

3751 

PI 

35 

0600 

31 

4159 

2751 

ST  ART 

R*U 

CH  K 

26 

1000 

60 

010  7 

0011 

BM  1 

BOOO 

27 

0011 

46 

6000 

0015 

R  A  U 

BONE 

38 

0015 

60 

010  3 

0007 

FSB 

1700 

39 

0007 

33 

1700 

0037 

N!E 

L0P1 

30 

0037 

45 

0030 

0031 

LDD 

EOOCL 

31 

0030 

69 

003  3 

0036 

PCH 

1977 

33 

003  3 

71 

197  7 

0077 

PCH 

1977 

33 

0077 

71 

197  7 

0137 

PCH 

1977 

L0P1 

34 

0127 

71 

197  7 

0031 

I  OP  1 

0  «  11 

R 

35 

0031 

60 

0105 

0009 

FMP 

8003 

36 

0009 

39 

800  3 

0013 

STU 

RSOD 

37 

0013 

31 

0018 

0031 

R  A  U 

NM  A  X 

3B 

0021 

60 

0100 

0005 

F  A  0 

ONE 

39 

0005 

33 

0050 

0177 

8TU 

DISKS 

40 

0177 

31 

003  3 

0035 

FIIP 

T  Z  RO 

4  1 

0035 

39 

010  1 

0001 

STU 

H  1  IE 

CYL  HEIGHT 

42 

OOO  1 

21 

0006 

0059 

FMP 

MU8 

4  3 

0059 

39 

0102 

0003 

STU 

CY  MFP 

C'V  L  MFP 

4  4 

0002 

31 

0056 

0109 

FAD 

BONE 

45 

0109 

33 

010  3 

0039 

8TU 

BTHRE 

4  6 

0029 

31 

00  3  4 

0037 

LflO 

ZR  0 

47 

0037 

69 

00  0  0 

0003 

STO 

N 

4  8 

0003 

24 

0156 

015  9 

STO 

AN8 

LOPS 

4  9 

0159 

34 

no  l  2 

0065 

LOPS 

R  AU 

N 

50 

00  65 

60 

016  6 

0061 

F  AO 

PT  4  5 

51 

0061 

33 

0350 

0237 

FMP 

TZRO 

53 

0227 

39 

ni  0  1 

0051 

STU 

0  1  8  T 

S3 

0051 

21 

020  6 

0209 

FMP 

MUS 

54 

0209 

39 

010  3 

0052 

F  AO 

BONE 

55 

0052 

32 

010  3 

0O79 

STU 

ARC 

56 

0079 

31 

00  8  4 

008  7 

FSB 

FP  100 

57 

0087 

33 

04  0  0 

0377 

BM  1 

BO  00 

58 

027  7 

46 

008  0 

6000 

R  AU 

AR  G 

59 

0080 

60 

008  4 

0039 

LOO 

EOHE 

60 

0039 

69 

004  3 

0045 

STU 

EONEX 

61 

0042 

31 

004  6 

004  9 

R  AU 

D  1  8T 

63 

0049 

60 

020  6 

0111 

FAD 

A 

63 

0111 

33 

010  4 

0081 

STU 

H 

64 

0081 

31 

0086 

0089 

FMP 

8003 

65 

0069 

39 

800  3 

0043 

STU 

HSOO 

66 

004  3 

31 

004  8 

0151 

R  AU 

R  800 

67 

0151 

60 

0018 

0033 

FOV 

H  S  OD 

68 

0023 

34 

004  8 

0098 

F  AO 

ONE 

69 

0098 

33 

OOSO 

0337 

LDD 

SOROT 

70 

0337 

69 

0130 

0083 

STU 

8E  CI 

71 

0130 

31 

013  4 

0137 

FMP 

ARC 

7  2 

0137 

39 

0084 

0184 

STU 

AROSE 

73 

0184 

31 

003  8 

0041 

FSB 

FP  100 

7  4 

00  4  1 

33 

04  0  0 

0377 

BM  1 

LOP  3 

75 

0377 

46 

0180 

0131 

R  A  U 

A  R  G  8E 

76 

0180 

60 

00  3  8 

0093 

LOO 

EONE 

77 

0093 

69 

009  6 

0045 

STU 

EONET 

L0P4 

78 

0096 

31 

065  0 

0053 

L0P3 

LOO 

ZR  0 

79 

0131 

69 

OOOO 

0153 

STO 

EO  NEY 

L0P4 

80 

0153 

24 

065  0 

0053 

L0P4 

R  A  II 

EONEX 

81 

0053 

60 

004  6 

0301 

F  8R 

EONEY 

83 

0201 

33 

06S0 

0437 

FAD 

A  N  S 

83 

0427 

33 

0013 

0139 

STU 

A  N  S 

84 

0139 

31 

001  3 

0115 

R  AU 

NM  A  X 

85 

0115 

60 

0100 

0055 

FSB 

N 

86 

0055 

33 

0156 

0133 

NZE 

LOPS 

87 

0133 

45 

0136 

0187 

R  III 

N 

88 

0136 

60 

0156 

0161 

F  AO 

ONE 

89 

0161 

33 

0050 

0477 

STU 

N 

L0P2 

90 

0477 

31 

0156 

0065 

LOPS 

RAU 

AN  8 

91 

0187 

60 

0012 

0017 

FMP 

S 

93 

0017 

39 

0106 

0356 

FOV 

TWO 

93 

0256 

34 

0150 

0700 

FMP 

TZRO 

94 

0700 

39 

0101 

0251 

STU 

FLUX 

95 

0251 

31 

0306 

0259 

STU 

1977 

96 

0259 

31 

197  7 

0230 

RAU 

H  1  TE 

97 

0230 

60 

000  6 

0211 

FMP 

P  1 

98 

0311 

39 

060  0 

0750 

FMP 

R  300 

99 

0750 

39 

0016 

0068 

STU 

VO  L 

CYL  VOLUME 

100 

0068 

21 

002  2 

0025 

FMP 

THREE 

101 

00  25 

39 

0300 

0800 

FOV 

P  1 

103 

0  800 

34 

060  0 

0850 

FOV 

RSOD 

103 

0650 

34 

0018 

0118 

FAD 

A 

104 

0118 

33 

010  4 

0181 

FMP 

A 

105 

0181 

39 

0104 

0004 

FMP 

A 

106 

0004 

39 

010  4 

0054 

LDO 

CBR  OT 

107 

0054 

69 

005  7 

0010 

FSB 

A 

108 

0057 

33 

0104 

0231 

STU 

H  P  R 

EOU 1 V   H 1 TE 

109 

0231 

31 

018  6 

0189 

FMP 

MUS 

110 

0  189 

39 

010  3 

0153 

STU 

MF  PPR 

111 

0153 

31 

0356 

0309 

F  AO 

BONE 

113 

0309 

33 

010  3 

0139 

STU 

|  3PR  M 

113 

0139 

31 

023  4 

0337 

RAU 

RSOD 

114 

02  37 

60 

00  1  8 

0073 

FOV 

A 

115 

0073 

34 

010  4 

0154 

-0  000  0  00 0  0000 00 0  0 
-000000000-000000 
-00 00 0  00 00-00 00 0  0 
-OOOO 0  00 0  0- OOOO 0  0 
-000000000-000000 
-000000000-000000 

-ooooooooo-oooooo 

-000000000-000000 

-ooooooooo-oooooo 
-ooooooooo-oooooo 

-000000000-000000 

-ooooooooo-oooooo 
-ooooooooo-oooooo 
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SOR  OT 
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SE  CL0 
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LOO 

ET«0 

STU 

PARI 

R  AU 

BONE 

FHP 

SE  CLO 

STU 

A  BC1 

F  SB 

FP  100 
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L0P6 

Rill 

A  B  C  1 
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SE  CLO 

S  TU 
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LOO 
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STD 
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L0P7 

L0P7 
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APPENDIX  B 

Description  of  IBM  650  Computer  Program 

Used  to  Calculate  the  Error  in  Centerline  Uncollided  Flux 

Determination  from  the  End  of  a  Finite  Right-circular  Cylinder 

by  the  Use  of  an  Equivalent  Circular  Plane  Source 

This  program  was  written  to  calculate  the  error  of  a  flux  determi- 
nation by  the  empirical  method  developed  in  this  thesis.  The  program  was 
written  in  Soap  II  computer  language  using  floating  point  operations.  The 
logic  diagram  and  the  object  program  are  given  in  this  appendix. 

The  solution  to  the  equivalent  circular  plane  source  approximation 
was  designated  ABC9  and  the  solutions  programmed  were 

*(a)  =  h-  {  Ejft);,)-^  (b2  sec  0z)j  ,   if  ^h  <  1.5;         (B-l) 

<p(a)  =  I  {  El(b2)  -  E1(b2  sec  6Z)}  ,  if  ^h>1.5.        (B-2) 

2M-a 

Equations  (B-l)  and  (B-2)  appear  as  Eqs.  A-6  and  A-7  in  Appendix  A.  The 
program  subroutine  for  E..  (x)  is  identical  to  that  described  in  Appendix  A. 
Input  data  and  output  data  associated  with  this  program  are  listed  in 
Tables  B-l  and  B-2  respectively. 

Table  B-l.  Input  data  required  for  use  of  the  IBM  650  Computer 
program  which  calculates  the  error  in  the  empirical 
solution  of  the  centerline  uncollided  flux  from  the 
end  of  a  finite  right-circular  cylinder. 

Symbol                     Explanation              Drum  Storage 
..  Location 

HITE  Cylinder  height  (cm)  0000 

BONE  Shield  mean  free  paths  0001 

FLUX  Actual  flux  (cm"2  sec-1)  0003 

S  Source  strength  (cm~3  sec-1)  0004 
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Table  B-l.  (concl) 


Symbol 


Explanation 


Drum  Storage 
Location 


R 

A 

MUS 

MUSZ 


Cylinder  radius  (cm)  0005 

Shield  thickness  (cm)  0006 

Source  cross  section  (cm--*-)  0008 

Empirically  determined  self-absorption  0009 
distance  (mean  free  paths) 


Table  B-2.  Identification  table  for  print  outs  for  error 
calculation. 


Card 


Word 


Symbol  and/or  Explanation 


1 
2 
3 
4 
5 
6 
7 
8 
1 
2 
3 
U 
5 
6 
7 
8 


FLUX  accurate  flux  solution 

ABC9  approximate  flux  solution 

ERROR    ABC9  -  FLUX 

PCT     (ABC9  -  FLUX)/FLUX 

z     self-absorption  distance 

A 

R 

BONE 

HITE 

MUS 

M^h  Cylinder  height  (mean  free  paths) 

p.  Shield  cross  section 

MUSZ 

MUSZ +  BONE 

A  +  R 

A/R 
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OBJECT  PROGRAM  FOR  CALCULATION  OF  ERROR  IN  EQUIV- 
ALENT CIRCULAR  PLANE  SOURCE  FLUX  DETERMINATION 
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An  IEM  650  Computer  code  was  written  to  solve  for  the  flux  from  the 
end  of  a  finite  circular  cylinder.  This  code  was  then  used  to  solve 
for  the  flux  from  cylinders  of  varying  size  and  absorption  coefficient  as 
well  as  for  varying  shield  thickness  and  absorption  coefficients. 

From  the  information  obtained,  a  method  was  developed  to  determine 
fluxes  from  cylinders  by  considering  them  to  be  replaced  by  a  circular 
plane  source.  This  circular  plane  source  is  located  within  the  confines 
of  the  cylinder  and  has  the  same  radius  as  the  cylinder.  This  method  makes 
it  possible  to  determine  the  uncollided  flux  from  the  end  of  a  cylinder 
more  accurately  than  has  previously  been  possible  with  only  a  desk 
calculator. 

Curves  and  supporting  information  as  necessary  are  presented  to  enable 
users  to  determine  the  source  strength  of  the  circular  plane  source  and 
its  location  within  the  cylinder. 


